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Resume. On demontre dans cet article un Theoreme de Lagrange, 
pour un certain algorithme de fraction continue en dimension 
2, dont la definition geometrique est tres naturelle. Des pro¬ 
prietes type Dirichlet sont aussi obtenues pour la convergence 
de cet algorithme. Ces proprietes proviennent de caracteristiques 
geometriques de l’algorithme. Les relations entre ces differentes 
proprietes sont etudiees. En lien avec l’algorithme presente, sont 
rapidement evoques les travaux de divers auteurs dans le domaine 
des fractions continues multidimensionnclles. 

Abstract. A Lagrange Theorem in dimension 2 is proved in this 
paper, for a particular two dimensional continued fraction algo¬ 
rithm, with a very natural geometrical definition. Dirichlet type 
properties for the convergence of this algorithm are also proved. 
These properties proceed from a geometrical quality of the algo¬ 
rithm. The links between all these properties are studied. In re¬ 
lation with this algorithm, some references are given to the works 
of various authors, in the domain of multidimensional continued 
fractions algorithms. 


1. Introduction and results 

1.1. Quick presentation of the main results. Since the beginning of 
the theory of Multidimensional Continued Fractions, an extension of the 
well known Lagrange Theorem in dimension one has been searched for. 
Historical remarks on the multidimensional continued fractions (the Jacobi- 
Perron algorithm and others) can be found in the works by F. Schweiger: 
[22] and [23], and A.J. Brentjes: [3]. 

The classical one-dimensional continued fraction algorithm applied on a 
real number x generates a sequence (£ s ),, gN in M, with £o = x i named the 
”complete quotients”, and Lagrange proved that the following assertions 
are equivalent: 

(1) x is a quadratic algebraic number. 
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(2) There exist natural numbers s > 0 and p > 1 such that £ s+p = f s . 

(3) There exist natural numbers so > 0 and p > 1 such that for every 
s > so, 6i+ P = £s holds (periodicity). 

The property (2) will be called loop property in this paper. 

Here we define a very natural two-dimensional continued-fraction algo¬ 
rithm for which the analogue in dimension two of the properties (1) and 

(2) are equivalent: This algorithm, named Smallest Vector Algorithm or 
” S VA” , makes a loop (property (2)) if and only if the real numbers which 
are its two initial values are in the same cubic field (property (1)). The 
SVA is defined at the beginning of Subsection 1.3.. 

We have to notice that we do not have periodicity, i.e. the property 

(3) , for initial values in the same cubic fields. The reason why is that our 
algorithm, unlike a lot of known multidimensional continued fraction algo¬ 
rithms, is not of the vectorial kind. Therefore, the loop property (2) does 
not imply periodicity (3). Nevertheless, the loop property (2) implies inter¬ 
esting algebraic properties and the fact that the algorithm is not vectorial 
permits strong approximation properties. 

Let’s state our Lagrange-type theorem. From any initial value 
Xo = X = r (xo, xi, X 2 ), with 0 < xo < x\ < X 2 , the Smallest Vector 
Algorithm generates a sequence (X s ) = ( r (xo, s , xi )S ,^ 2 ,s)) of triplets of 
real numbers, and we have the following statement. 

Theorem 1 (Lagrange Loop Theorem). First Part: Let p be any real 
root of a third degree irreducible polynomial P (r) = r 3 — ar 2 — br — c, 
with a,b,c rationals; let X= r (xo,xi,X 2 ) be any rationally independent 
triplet of real numbers in the field Q [p], with 0 < xq < x\ < X 2 - Then the 
Smallest Vector Algorithm applied on the triplet X ”makes a loop”: there 
exist integers s and p with p > 0 and a real number X such that: 

X s+P = AX S or equivalently: x s+p = x s , with x s = ( ——, —^ 

\X2,s x 2 , s 

Moreover, A is an algebraic integer of degree 3, and a unit, such that 
Q [p] = Q [A]. The minimal polynomial of A can be easily deduced from 
the relation X s+p = AX S , as also the expressions of and as rational 
fractions of X. 

Second Part: Converse Statement: Let X = r (xo,xi, X 2 ) be any ra¬ 
tionally independent triplet of real numbers, with 0 < xo < x\ < X 2 - Let's 
suppose that the Smallest Vector Algorithm applied on the triplet X makes 
”a loop” i.e. that X s+p = AX s with p > 0. Then A is an algebraic integer 
of degree 3, and a unit. Again, the minimal polynomial of A can be easily 
deduced from the relation X s+p = AX S , as also the expressions of and 
as rational fractions of X. 
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The objects in this theorem are more precisely described in the follow¬ 
ing subsections. We also prove that the same algorithm provides rational 
approximations with Dirichlet properties, id est, with an optimal exponent. 

Throughout this paper, we are going to use only the canonical euclidean 
norm and inner product in M 3 for our approximations. 

Our Dirichlet property is that for every independent triplet of posi¬ 
tive real numbers X = 1 (xq,xi,X 2 ), the algorithm generates a sequence 

^g ( < ) s l . g 1 ^, g|> s ^ of triplets of three-dimensional integer vectors, which real¬ 
ize integer approximation of the plane (X 3 -) with the following inequality, 
on an infinite set S of integers: 


sup 

seS 


( min 

g| s) • X 

\z=0,l,2 1 



( max 

y i=0,l.,2 



2 


< Too 


(the index ^ is above, in parentheses; the big point denotes the scalar 
product), with of course: lim ( max ||g. s ^||^ = Too. See Theorem 


2 in subsection 1.3.. We prove additional Dirichlet properties, for the in¬ 
teger approximation of MX as well as of X^, when X^ (or X) has a bad 
approximation property. See again subsection 1.3.. 

Let’s notice that the approximation properties of our algorithm hold only 
for a subsequence of the integer vectors (gi A ' ) , g|' s) , g^ ) ; the algorithm 

has a very simple geometrical definition, and strong geometrical, algebraic 
and approximation properties, but it is not designed to provide only best 
approximants, or only approximants with optimal exponent. 

The goal of this paper is also to show the relations between different 
kinds of properties of such an algorithm: 

(a) Lagrange property, (b) Dirichlet approximation properties', 

(c) Best approximation properties. 

(d) Properties of the triplet X which is the initial value of the algorithm 
(it may be badly approximate by integers, or well approximate) 

(e) Geometrical properties of the tetrahedrons formed by the three integer 
vectors, generated at each step by the algorithm. 

This study is a generalization of the well known continued fractions the¬ 
ory in dimension 1, with a two-dimensional algorithm which has more prop¬ 
erties than most of the existing ones. Let’s notice that all the mathematical 
techniques used in this paper are elementary. The most sophisticated tool 
appearing here is the Minkowski’s Theorem on Successive Minima of sym¬ 
metrical convex sets. 

At the end of the paper, in Section 6., the reader shall find a short 
review of the themes on Diophantine approximation which are related to 
this paper, with some bibliographical references. 
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In subsection 1.3. are given the main theorems and definitions of this 
paper. In subsection 1.4. the reader shall find two numerical examples of 
Lagrange loops. The plan of our paper is in subsection 1.5. 

But first, in order to understand better the two-dimensional case, we 
recall some facts and notations about the classical continued fractions al¬ 
gorithm in dimension one, from a particular point of view. 

1.2. The one dimensional example. 


1 . 2 . 1 . A formalism with matrices. A real number x is chosen, with 0 < x < 
1. Here it is supposed to be irrational. Let the vector X be: X = T 0,l), 
where the ” T ” denotes the transposition. The classic one-dimensional 
algorithm provides integer points T (p n -, ( ln) which are the nearest inte¬ 
ger points to the line ID) = MX. These points are called ’’convergent” 

is _ fPn—i Pn\ yy e h ave the re- 


points. We consider the matrices B n = 


lation: B„+i = B, 


, where a n is the ra-th ’’partial quotient” 


of the continued fraction and is a strictly positive integer. If we denote: 

A n = ( ), then the approximating matrices B„ appear as products 

of matrices A*. (1 < k < n— 1). Let’s notice that Bi is the Identity matrix. 

In order to be closer to our two-dimensional algorithm, we may also 
split the n-th step into more elementary steps, and consider the simple 


matrix D = 


then we have B„+i = B„A„ = B„D a 


The last matrix corresponds to an exchange of vectors, when the following 
convergent T (p n+ i,q n+ i) is found. We may notice that all the matrices 
involved have determinant ±1. 


1.2.2. Polar matrices, cofactors, periodicity. We also introduce the polar 
matrices G n , each of them being the transposed matrix of the inverse of B„ . 
Let g 0 ,„ and gg n be the column vectors of G n . These vectors realize integer 
approximations of the line A orthogonal to B, and the regular continued 
fraction algorithm is precisely designed to obtain both: g 0)n • X > 0 and 
gi )T( , • X > 0 for the scalar products. 

These quantities gy n • X and gi jn • X are particularly important in the 
theory of continued fractions. Let bo. n and bi )Tl be the column vectors of 
B n . We have the vectorial relation: (g 0>n • X) b 0;n -|- (gi, n • X) bi jn = X. 
(To see that, make the scalar product of the left-hand vector of the equality 
with go,n> and then with gi jn ). Because of this relation, (go,n • X) and 
(gi >n • X) are called the cofactors in the algorithm. 


Let’s form the cofactors vector: X n 


f go,n • X \ 
\Sl,n»Xj 


. We have: 
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r G n X = X n , id est B n 1 X = X n . Now we may calculate X n . 

We have B” 1 = (—l)( n_1 l f ^ n and then 

n \~Qn-l Pn-lJ 

x n = T ((-l) n (Pn - xq n ), (-l) (n-1) (p n _ 1 _ xq n _ i)^. 

Then we define the sequence (x n ) by x n = (—l) n (p n — xq n ) = go, n • X. 
At the first step, x\ = x. 

We may now give the rule which defines the quantity a n , in dimension 


1. Here the brackets denote the integer part: a n = 


gi. 


X 


Lgcu 


X 


%n— 1 


Xr, 


We now define another object, the inverse of this quotient: := 

Then the preceding relation writes: a n = 

The quantities x n and the cofactors vectors X n are of highest interest 
in questions concerning Lagrange property. This algorithm is eventually 
periodic, from the range n, if and only if there exist an integer p > 1 and a 
real number A such that X n+p = AX n , or equivalently: £ n +p = £n- These 
are the conditions we shall use in our Lagrange Theorem. 


1.2.3. Recursive relations on the polar matrices. Let’s denote by M* the 
polar matrix of M, such that M*=( 1 M) 1 . We have G n = B*, and 
then, each matrix B n is a product of matrices A£ (1 < k < n — 1), with 

a; = (T Ij . We may consider the simpler matrix C = ^ 


1 0 

\ / 

then we have G n+ i = G n A* = G n C 


-1 1 


0 1 
1 0 


This leads to the relations 1 G n+ \ = A^ 1 1 G n and then, by T G n X = 

Xn, to. X„. | — A n X n — ^ j Xn, therefore. x n +i — x n —i a n x n . 

Because of this formula, the continued fractions algorithms may be called 

i 


subtractive. This leads to the recursive relation: £ n +i = -J- — a n = J- — 

S n 

In particular, if go and gi are the column vectors of G, then the co 
vectors of the following matrix G x C = G x (\ !) are (g 0 - gi) and 


f n 

umn 


. b 

gi. Our two-dimensional algorithm is built in a similar way, in the next 
subsection. 


1.2.4. Use of the orthogonal projections. From a more geometrical point 
of view, let’s denote by gg n and g" n the orthogonal projections of go, n and 


gi >n on Hi. Then we also have: X n = ||X| 


60 ,r 
g l,r 


Concerning a Dirichlet property of the algorithm, it can be written: 
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For each n, ||gi „|| Qn A 1- Our matricial formalism would permit to 
prove easily the Lagrange Theorem for a quadratic number t, using this 
Dirichlet property. 

Now we are going to generalize all these properties and demonstrations 
in dimension two. 


1.3. Results in this paper. Now we are in dimension two. Throughout 
this paper we suppose that the triplet X = r (.To, xi, £ 2 ) of real numbers is 
rationally independent, and verifies 0 < To < Ti < £ 2 - We use the canonical 
euclidean norm in R 3 and the canonical scalar product. We denote by ID 
the line ID = lx and by P the plane orthogonal to ID. 

As it is usually done in this field, the index ^ of the sequences will be 
above, in parentheses. 


Definition. The Smallest Vector Algorithm is described by the following 
sequence (G^ s l) sgN of 3 x 3 integer matrices, which is inductively defined 
by: 

a) G® = I, the Identity matrix. 

b) Let’s suppose G^ = ^gjj s \ g^, g^ ^ has been defined. Let g,-^ and 
denote the respective orthogonal projections of gf 1 on D and P. Let 


A min denote A min := min ( gf s) - gy' s) 


l(s) t(s) 

g 2 - Si 


l(s) l(s) 

g 2 - So 


We define first the three column vectors (fo,fi,f 2 ) of G^ s+1 \ in disorder: 


(fo,fi,f 2 ) 

= (go S) . 

(s ) 

gl 

(s) 

-So 

(s)' 

■S 2 

j if A m i n = 

r(s) 

gl - 

l(s) 

So 

1 

(fo,fi,f 2 ) 

= (go S) . 

(s ) 

gl ■ 

(s) 

g 2 - 

(sy 

gl 

j if A m i n = 

t(s) 

g 2 - 

/(s) 

Si 

1 

(f 0 ,fl,f 2 ) 

= (so\ 

( s) 

gl ■ 

(s) 

g 2 - 

(s)' 

go 

j if ^min = 

t(s) 

g 2 - 

/(s) 

go 


g( s+1 ) = 

U s+1) 

(s+l) (s- 

. gl . g 2 

«,)' 

is defined ae 

1 any of the 

permutations 

of the vectors (Fq , f 

1 , h) 

such 

that 

we have: 

//(s+l) 

g 0 

< 

g'f +1) < 


//(s+l) 

& 2 


The columns gg S \g^,g^ of the matrices G^ realize integer approx¬ 
imations of the plane P. They play the same role in dimension 2 as, in 
dimension 1, the matrices G n studied above. 

As in dimension 1, the cofactors vector X^ is fundamental. It’s defined 
by: X( s ) = ||X|| • 1 ^ g^ 6 ' 1 , g 7 /^ , g^^ ^ , and the vector the 

/ „"<•) „"(•) \ 

’’projective” version of X (s ), is defined by = 1 ■■ ) ■■, ■■ . ■■, 1 . 

y s 2 s 2 J 

The Smallest Vector Algorithm has a Lagrange property, which is ex¬ 
pressed by Theorem 1 of subsection 1.1. The demonstration of this theorem 
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(see below) is essentially the same as in dimension one, based on a Dirichlet 
result (In dimension two, Theorem 2). 

First we introduce the unimodular positive integer matrix defined 
as the polar matrix of (the transposed matrix of the inverse of G^). 
We denote by bg , b) , bg the column vectors of which realize integer 
approximations of the line ID. The vectors b^ 5 ) will be the 

orthogonal projections of these vectors on B, and the b ; ' their orthogonal 
projections on P. In the same way are defined the gjh 5 ) and the g} s \ and, 
for any vector h, the projections Id and h" of h on P and B, 


Theorem 2 (Dirichlet Properties). For each rationally independent triplet 
X of real numbers, with 0 < xo < x\ < X 2 , the Smallest Vector Algorithm 
has the following properties 

a) There exists an infinite set S of natural integers such that: 


sup 

sGS 


( max 
\i=0,l,2 



2 


( min 

V=0,l,2 



< Too, with also: 


lim ( min llgf^ll 1=0 .In other words , this alqorithm provides a 

s->+oo, seS \i=0,l,2 11 1 "y 

Diophantine approximation of P which possesses a Dirichlet property, with 
the optimal exponent, 2, and with a form which is rather strong, 

b) IF there exists c > 0 such that for any integer point h / 0, 


|h|| 2 ||h"|| > c holds (id est, if the couple 


xi x 2 
x 0 ’ x 0 


is badly approximable, 


as it is proved in Section 3), then the following relation holds with the 
Smallest Vector Algorithm, with the same infinite set S: 

r / 2 / \ i 


sup 

seS 


[ max 

g ; (s) 

V i=0,1,2 



max 

*= 0 , 1,2 



< Too, with also: 


lim ( max ) =0. 

s-H-oo, seS \*=o,i,2 11 M y 

c) If again there exists c > 0 such that for any integer point h / 0, 
||h|| 2 ||h"|| > c holds, then 


sup 

ses 


( max 

\i=0,l,2 



( max 

b" (s) 

)1 

\i= 0,1,2 


)\ 


< Too; with also: 


lim 

( 

max 

bf } 

s — s^S 

1 *=0,1,2 



Theorem 3 (Geometrical Theorem). The Smallest Vector Algorithm pos¬ 
sesses the following property. Let be the convex hull 
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H'W = 


CONV 


/(«) '(*) _'(s) 


'(a) 
-gl ■ 


'(a) 


}) * nP - Let be 


(( '(a) '(a) /(a) /(s) 

(|g 0 ,gl >§2 >-S0 .“gl ~.g2 
f/ie radius of the greatest disk in P with center 0 contained in H h s l. Then 
there exists an infinite set S of natural integers such that: 

'(a) I 


sup- 

ses 


max 

#= 0 , 1,2 


g# 


P 


(a) 


< +oo. 


max; g. 

Definition. If lim inf l_0 ’ 1,2 


/(»)! 


o(«) 


< +oo, with the notations of the above 


s —^-(-oo P 

Geometrical Theorem, it will be said that the algorithm is balanced. 


As a Best Approximation result, we give our Prism Lemma, which is 
proved in Subsection 3.1.. It is a very easy result, but it is true and 


important for any continued fraction algorithm, and the author has not 


seen this statement anywhere in literature. 

Lemma ( Prism Lemma ). Let g^\ be the column vectors of 

the matrix Gl s ) generated at the s-th step by the Smallest Vector Algorithm. 


Let the sets be the convex hulls: 

H /(s ) = CONV 


'(a) 

§0 . 


'(a) 

gl • 


/(a) 

g 2 > 


/(a) 

-go > 


/(a) 

-gl > 


'(a) 

-go 


in P . 


We shall omit the indices ^ s \ Let H be the prism: H = H' + B. Then, 
with the usual notation h" for the orthogonal projection of h on B: 

-For each non zero integer point h in H , ||h ,/ || > llgoll holds. 

-For each integer point h in H which is not of the form nogo, with no 
integer, ||h ,, || > 11g"11 holds. 

-For each integer point h in H which is not of the form nogo + nigi, 

with no and n\ integers, ||h ,/ || > Hg^H holds. 

- That implies that, if (ho,hi,h 2 ) is a free triplet of integer points in 

H, then max (||h''||) > max (||g"||) = ||g 2 1|- 

i=0,1,2 V *=0,1,2 VM 


The previous theorems suppose that the initial values (xq\x\-,X 2 ) of our 
algorithm are rationally independent. The reader may wonder what hap¬ 
pens when they are not. 


Theorem 4. The Smallest Vector Algorithm finds rational dependence. If 
the initial values (xq,xi,X 2 ) of our algorithm are rationally dependent, then 


the SVA generates at some step an integer triplet 


(. s) (s) (s) 

go >gl >g2 


such that: 


//(a) 

go 


.(*) 


= g 0 ' • X = 0. This relation gives the coefficients of the rational 
(integral) dependence of (xq, £i, £ 2 )- 


The proof uses Lemma 2 in the next subsection, the Geometrical Theo¬ 
rem and the Prism Lemma. The demonstrations of these three results are 
still valid if (xo,x\,X 2 ) is rationally dependent. With these three results, 
the proof of Theorem 4 is very short. 
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Proof. Let’s suppose that we have an integral dependence relation, of the 
shape h • X = 0, h being a non null integer vector. By the Lemma 


( max 

/(s) 

g i 

W=0,l,2 



2 of the next subsection: lim 


by the Geometrical Theorem above, lim inf 

s->+oo 


max 
i= 0,1,2 


= Too. 
.'<•>11 


In addition, 


Si 


o( s ) 


< Too. Then 


lim sup ff s ' > = Too. Then h belongs to some hexagon H'W defined as 

s—>+oo 

a convex hull in the Prism Lemma just above. By this Prism Lemma, 


0 = ||h"|| > ||gg|| = 


x. 


So 


(s) 


IXII 


> 0 holds. Then X 


So } = 0 . 


□ 


1.4. Numerical examples. We give two examples of Lagrange Loops 
when the initial values are in some cubic field. 

Example. (xo;xi;x 2 ) = (l;2cos (ir/ 7) ;4cos 2 (7r/7)). The sequence of the 
is the simplest the author has met. Almost every couple in it 


M (s) 
Xg m x\ 

r (s) ’ O) 


repeats infinitely. We have: xf — x\ — 2x\ + 1 = 0, and 

(x ‘2 — 2) 3 T (x 2 — 2) 2 — 2 (x 2 — 2) — 1 = 0; the algorithm provides the fol¬ 


lowing ( s T 1, 


(s) (s) 

x 0 . X 1 
r ( s ) ’ _(«) 
•"2 2 


{1, {1-80193773580484, 
{2, {1.80193773580484, 
{3, {1.24697960371747, 
{4, {1.24697960371747, 
{5, {1.80193773580484, 
.../... 

{292, {1.24697960371747, 
{293, {1.80193773580484, 
{294, {1.24697960371747, 
{295, {1.24697960371747, 
{296, {1.80193773580484, 
{297, {1.24697960371747, 
{298, {1.80193773580484, 
{299, {1.24697960371747, 
{300, {1.24697960371747, 


3.24697960371747} } 
2.24697960371747} } 
2.80193773580484} } 
1.80193773580484} } 
2.24697960371747} } 

1.80193773580484}} 

2.24697960371747}} 

1.80193773580484}} 

1.55495813208737}} 

2.24697960371747}} 

1.80193773580484}} 

2.24697960371747}} 

2.80193773580484}} 

1.80193773580484}} 


Example. (xo]X\;x 2 ) = ( 1; v^ 13; 13 2 ). The algorithm provides the fol 


Here there are very few repetitions 


lowing sequence T 1 
(loops), but they exist. 


«(») T W 

(s) 1 (s) 

x 2 X 2 


{1, {2.35133468772075748950001, 5.5287748136788721414723}} 
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{2, {2.35133468772075748950001, 
{3, {1.35133468772075748950001, 
.../... 

{104, {2.35133468772075748950001, 
{105, {2.35133468772075748950001, 
{106, {1.35133468772075748950001, 
../... 

{160, {5.5057084068852398563646, 

.../... 

{219, {1.35133468772075748950001, 

.../... 

{308, {2.35133468772075748950001, 

.../... 

{411, {2.35133468772075748950001, 


4.5287748136788721414723}} 

4.5287748136788721414723}} 

5.528774813678872141472}} 

4.528774813678872141472}} 

4.528774813678872141472}} 

47.82563987973201144317}} 

4.5287748136788721414723}} 

5.5287748136788721414723}} 

5.5287748136788721414723}} 


1.5. PLAN of the paper. We shall prove the Results above in the order 
in which they are written in this first section, and, in fact, in the reverse of 
the logical order. 

• In the second section, we admit that the Smallest Vector Algorithm 
or SVA has the Dirichlet Properties of Theorem 2, and we prove 
that this implies the Lagrange properties. But this demonstration 

is made only in a particular case, namely X = ^1, \ZlV, \/N 2 ^j 
with N natural. The complete proof of the Lagrange Property for 
three numbers in a cubic number field is given in Section 5. But 
this general demonstration is intricate, and the particular case gives 
all the main ideas involved. That’s why we prefer to begin with this 
particular case, for more clarity. 

• In the third section, we admit the geometrical property of the SVA, 
namely that it is balanced. We prove that this property implies the 
Dirichlet Properties. 

• In the fourth section, we prove the Geometrical Theorem, namely 
that the SVA is balanced. 

• In the fifth section, we give the general demonstration of the La¬ 
grange Theorem. Then all theorems are proved. 

• In Section 6, we locate the results of this paper in relation to the 
main themes in Multidimensional Dimensional Fractions Theory and 
Homogeneous Diophantine Approximation, with some bibliographi¬ 
cal references. 


2. Demonstration of Lagrange Properties (particular case) 
2.1. Four basic Lemmas. We’re going to need the following lemmas. 
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Lemma 1 (Basic Properties of Brentjes’ Algorithms), a) Inductively, by 
the way of building the Smallest Vector Algorithm, we have the following 


properties: 

0 < x.g[, s) = ||X||||g''(*)|| < x.gj s) = ||X||||g'/W|| 
l|X||||ggW||. 

b) The following equality holds for every s € N: 


< X 


( s) 

1 §2 = 


IXI 


go ,(s) 


bl s) + IIXI 


gi' (s) 


b. s) + IIXI 


g' 2 /(s) 



= X 


That’s the reason why ||X|| ||go ( ' s - ) || , ||X|| ||gi ( ' S - ) || , ||X|| ||g"^|| are called 
cofactors. 

c) For every s € N: T G^X = X^ s ) or, which is the same: 

(B( S )) _1 X = X^). 


Proof. First of all, the a) Property is true at the step s = 0. Let’s suppose 
it’s true at the step s. The construction of the new gj s+1 ^ by subtraction 
is always done in accordance with the order of the X • g|) s ' i = ||X|| |g"(T||. 
Then, at the next step, each of the X • g^ +1 ' ) is positive, and then equals 
||X|| ||g( / ( s + 1 ) ||. The rule of the algorithm is to order these numbers at step 
(s + 1). Then the property is obtained at step (s + 1), and a) is true by 
induction. For the b) property, the equality 

^X • gQ S ^ + ^X • g^ •b^ + ^X • g^) -b^ = X holds. To see that, 
make the scalar product of both the left-hand and the right-hand vector of 

(s) 

the equality with each of the g) . 

This equality is the same as the equality in b). 

Because X^) = 1 ^X • gg S \ X • g^, X • g^), the equalities of c) are 
obvious. □ 


Lemma 2. For the Smallest Vector Algorithm, and more generally in any 
Brentjes’ algorithm, we have: lim ( max llgf'll I = Too. 


Proof. If this limit does not hold, there exist a real number M and 


infinite set T such that for any s € T, and for i = 0,1,2, 


g 


.'(«) 


an 

< 


M holds. In addition, by the subtractive nature of these algorithms, 
the ||g''( s )||, i = 0,1,2 are also bounded. Then the set of the triplets 

fg()^ ■ g] S ' 1 ■ gj) 5 ' 1 ') with s € T is bounded. Then it is finite. Then there 


exist two distinct natural numbers s and t such that ^g^, g^, g^' j = 

(gi'tgf.g?) and particularly: (gJM,g'/M, g?<*>) = (gg<‘>,gf«>, *£<*>). 
But that’s impossible, again by the subtractive nature of these algorithms. 
Our hypothesis was false, and then the conclusion of the theorem holds. □ 
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Now, we state two lemmas which will provide the Second Part of our 
Lagrange Theorem. 

Lemma 3 (Degree in a Loop). Let X= T (to, x\, t 2 ) be any rationally 
independent triplet of real numbers, with 0 < xq < x\ < t 2 . Let’s suppose 
that the Smallest Vector Algorithm applied on the triplet X ’’makes a loop ” 
i.e. that X^Lp) = AX(P with p > 0. Then A cannot be a quadratic real 
number. 

Proof. If the case of such a loop, we have X^ s+: p) = Axi s \ or equivalently 
(BP+p ))" 1 X = A (B ^)) -1 X or (B( s+ p )) _1 B^X^ = AXW. 

Let’s denote: B := (bP+p ))" 1 bW. Then: B X< s ) = AX^, with 

x (s) = T (4 s) ,xj s) ,4 s) )- 

Let’s denote: Y : = T := T (yo, 2 / 1 , 1 )- 

\ x 2 x 2 J 

Then also: BY = AY, which can be written 
/ Poo P01 P 02 \ / yo \ / Ay 0 \ 

Pm P 11 P12 U\ = Aj/i , where each /% is an integer. 

\ P 20 P 2 I P 22 / \ 1 / \ A J 

Then we have: fooVo + P 21 IJ 1 + P22 = A. From now on, we suppose that 
A is a quadratic real number. 

First Case: A 21 = 0 . In this case, yo belongs to Q(A). In addition, 
the same equality of matrices provides Pioyo + Piiyi + P 12 = Ayr, which 

can be written: Pi 0 yo + (fin — A) yi + P 12 = 0. Then yi = ^ 10 ^° + ^ 12 an d 

A — PH 

yi also belongs to Q (A). 

Second Case: p 2 i A 0. Then: y\ = 70 yo + 7 iA + 72 , with 70 , 71 , 72 ra¬ 
tional. The same equality of matrices provides Pooyo+Poiyi+Po 2 = Ayo, or 

(Poo + Poi 7 o - A)y 0 + /3 0 i7iA + 72 + A)2 = 0 . Then y 0 = l3 ° ni ^ + 72 + ^° 2 , 

A - Poo - PoiTo 

Then, in both cases, both yo and yi belong to Q(A), which is a vectorial 
space with dimension 2 over Q. Then the three numbers yo, yi and 1 are 
rationally dependent. So are therefore the three numbers Xq S \x[ s \x 2 S \ 
and then the three numbers xo = Xq°\ x\ = x[ 0> , t 2 = x^\ because 
X ( °) = B^X^P This contradicts our hypothesis. Then A cannot be a 
quadratic real number. □ 

Lemma 4. (Converse Statement in Lagrange Theorem). 

Let X = 1 (.To, x\, X 2 ) be any rationally independent triplet of real numbers, 
with 0 < To < ti < t 2 . Let’s suppose that the Smallest Vector Algorithm 
applied on the triplet X makes a ’’loop” i.e. that X^Lp) = AX( S ) with 
p > 0. Then A is an algebraic integer of degree 3, and a unit. The minimal 
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polynomial of X can be easily deduced from the relation = XX.^ as 

also the expressions of an d ix as rational fractions of X. 

Proof. By hypothesis,we have: X^' + p) = AX( S ), or equivalently 

(B^+p))" 1 X = A (B^) -1 X or (B (s+ p)) _1 (bW) X^ s ) = AX^. 

Let’s denote: B := (B^+p))” 1 (B^). Then: B X^ = AX^. Let F (£) 
be the polynomial: F (£) = det ^B — £1^, we have: F (A) = 0. But B is an 
integer matrix with determinant ±1; then we have a relation of the shape: 
A 3 + m\ 2 + nX ± 1 = 0, with m and n natural integers. 

Then A is an algebraic integer and a unit. Then, if A is a rational 
number, it has to be: A = ±1. But that’s impossible because we should have 
X( s+ p) = ±X( S \ which is impossible by the subtractive form of the recursive 
relation on the sequence (X^). By the previous Lemma A cannot be a 
root of a second degree polynomial. Then F (£) is the minimal polynomial 
of A over Q, and A is a cubic algebraic integer. Its norm is 1, and the 
relation A (A 2 + mX + n) = =Fl shows that A is a unit. 

Now, we have to solve the equation ("b — Al) X^ = 0, or (b - ai) Y = 


0, Y being the unknown vector, with the classic method of linear algebra. 
Because Y = (B^) X, the triplet Y is rationally independent, then 
2/3 / 0. If we just want to obtain one vector solution, we may even suppose 
that 2/3 = 1- Let A be the matrix ^B — Alj without its third row, let 

a 2 be the third column of the matrix A, and let A 7 be the matrix A 
without its third column. Let also Y 7 be the vector Y without its third 
coordinate. With block submatrices, we have to solve in Y 7 the equation: 

TY 7 1 
1 

J 

a 2 . We know that det (A 7 ) / 0, because F (£) is the 

jji-a (A 7 '" 1 


A 7 ;a 2 


= 0. This gives A 7 x Y 7 + a 2 x 1 = 0, or 
x J 

Y 7 = - (A 7 

minimal polynomial of A. Then we have X (s ) = 
a, and then: 


»2 


a 


, for some 


X = B (S )XW = qB< s ) x 
a multiplicative coefficient, 
the unit A. Then the vector Z : = 


-(A 7 )" 1 


a 2 


1 


Of course, X is defined up to 


Vote that A 7 and a 2 are rational fractions of 
- (A 7 ' 1-1 


a 2 


1 


can be also obtained as an 


expression of A. Let bo, bi, b 2 be the three rows of the matrix B^l; then 


we have 



Z = 



Then si = feoZ and = biz 

x 2 b 2 Z x 2 b 2 Z 


and we can express and as rational fractions of A. 


□ 



14 


Christian Drouin 


2.2. Demonstration of the Lagrange Property (particular case). 

In this subsection, we prove the Theorem on the Lagrange Property, admit¬ 
ting the Dirichlet Theorem which is proved in other section; we’re doing 
this proof only in the special case X = 1 (l,0,0 2 ), with 9 = \/~N, 

N being an natural number, but not 9. For the complete proof, see 
Section 5.. We’re going to prove that the Smallest Vector Algorithm in 
this case makes a loop. 

It is a basic result in Diophantine Approximation that there exists a real 
number e > 0 such that for any non-null integer point h = T (m, n,p ), the 
inequality (max (\m\ , |n|, |p|))" x |li • X| > e holds. See for instance [ 6 ] 
(Cassels), Theorem III, page 79, statement (2), which is much stronger. But 
in finite dimension, all the norms are equivalent. Then there exists d > 0 
such that for any non-null integer point h, the inequality ||h|| |h • X| > d 
holds, id est ||h||" ||h"|| x ||X|| > d. 

Then, by our Dirichlet Property Theorem 2, which we admit temporarily, 
there exists an infinite set of integers S such that: 


sup 

ses 


( max 

b' (s) 

Y( 

max 

b"W 

)1 

\2=0,1,2 


J \z=0,l,2 

z 

j\ 


lim ( max ||b'( s )|| ) = 0 . 

s— H-oo, sGS \i=0,l,2 


= L < Too, with in addition: 


Let s be any element of S, and let , b, , bl, j be the column vec¬ 
tors of (B*”)). We choose one of those three vectors, say b [ ) s ' 1 , which 
will be more simply denoted: b( s ) = b ( ' ) s ' i . Let’s define its coordinates: 
b*”) = r b{f\ bi s ^ . From now on and for a while, we may omit the 
indices . 


Notation 1. The notation M^ * = Mb or M [b^] will denote the integer 


j b z by b x \ 

matrix: Mb = M [b( s l] = \ N ■ b x b z b y , which has the inter- 

\N -by N ■ b x b z J 

esting property: MbX = Mb Id ] = (£u T b y 9 T b x 9 2 ) I 9 ]. Let’s 

W ' ' \ 0 2 ) 

denote by Ab or A [b^^] the following element of Z [0]: 

Ab := [b z + by9 + b x 9 2 ) ■ Then we have: MbX = AbX; i.e. Ab is an 
eigenvalue of Mb with eigenvector X. 


We consider the sequence of the matrices 11*”' = II = 1 MbG. Let 
(b##, b#, b) be the three column vectors of Mb- 
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b## • g 0 b#+ • gi b#+ • g2 

Then II = T M b G = { b* • g 0 h* • gi b # • g 2 J := ( 7 r,j), i = 

b • go b • gi b • gi 

0 ,1,2; j = 0,1,2. 

Lemma 5 (Main Lemma for Lagrange). The matrices 11+) are bounded 
independently from s. 

Proof. We have to find an upper bound for each of the | 7 Tjj|, i = 0,1,2; 
j = 0,1, 2. Let’s consider first the |b+ • g,; | . 

b x , 

= Q#b, with Q# = 


We have: h# = Q# 


X 



We also 


have Q+X = 0X. Let n be: n =-—-. Then also Q+n =0n. 

IIX.II 

We have, with always the same kind of notations, b+ = b #/ + b +". 
b # • g, = (b #/ + b # ") • (g' + g") = b #/ • g' + h*" • g", 

with g' = ± (b" A b'. + b' A b"); g" = ±b' A b'. . Then: 


b # *g, < b #/ b "j\ 

IK ||+ b#/ 

IK IIII 

b "k 

+ b # " 

ikiiimi- 

Let’s denote: B' = max 
*= 0 , 1,2 

b' | and P" = 

max 
2 = 0 , 1,2 

b" 

, so that: 

(ftf 13" < L 


Then: 

( 2 . 1 ) 


b # 


< 


b*' 


P"P' + 


b*' 


P'P" + 


b*" 


K) 


,l\2 


But we also have: b+ = Q+b = Q+(b^ + ||bo||n) = Q+bg + ||bo||0n. 
This proves first that the distance 11 b^ 7 11 between b+ and B is less than 

||Q #b o|| : 


( 2 . 2 ) 


b #' 


< 


Q#b' 


< 


Q 


# 


x INI ^ 


Q 


# 


P’ 


(this using the norm of the matrix). Furthermore, we have the following 
equality: b#' + b+" = b# = Q+bg + ||bg|| On, and we deduce: 
b#" = 11by 11 0n+ Q#b^ - b # '. Then: 


(2.3) 


b*" 


< bo 6 + 2 


Q 


# 


|b()|| <p"e + 2 Q* p ' 


Putting 12.21 and 12.31 in l2Tl we obtain: 


t)#(s) • g 
and then: 


W 


< (/?'+)) W (s) (2||Q # || +0) +2||Q#|| (0'W) 


b#M • g J s) < L (2 ||Q#|| + 0) + 2 ||Q#|| (/ 3 '+)) a . 


The limit of the last term is 0, by the c) of the Dirichlet Properties 


Theorem. Then the set of the 


b#+) 


g 


(s) 


, with s in S', is bounded. By a 
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similar demonstration, we prove that the numbers b^^ s ) • are also 

bounded and so are, in an obvious way, the numbers bg • g' : . Then the 

set of the 11^^ is bounded, and the Lemma is proved. □ 

With this Lemma, the demonstration of the Lagrange result is easy. 

Proof of the Lagrange Result. The set of the Il( s ) with s in S is bounded; 
but these matrices have integral coefficients. Then the number of all the 
nM is finite. Then there exist s and t, with s < t, such that 11^ = 11^). 
This means: i M [b^] x = 1 M [b^] x G^. By transposition: 

(B^) 1 x M [b^] = (B( s )) 1 x M [b^]. We apply that to the column 

vector X; we obtain (B®) 1 x M [bW] x X = (B ^) -1 x M [bW] x X; 
then, by ’’eigenvalue”, see Notation above: 

(B (t )) _1 x (A [bW] • X) = (B (s )) _1 x (A [bW] • X). 

We recall that A [b^ s ^] is a real number in Z [0]. 

Then (B ^) -1 X =^-[^ 7 TT ■ (B ls) ) _ 1 X. 

V ’ A [b«] V ’ 

This reads X (t ^ = XXf s \ with A = , and the first part of the 

Theorem is obtained. 

We have already proved the Second Part of the Theorem, by Lemma 4 
of the previous subsection. With this Second Part, we also obtain the final 
assertions of the First Part. 

In particular, xq,xi,X 2 , are rationally independent, then so are also , 
2i and 1. Then , ll is a basis of Q (p) over Q. But we know by the 

Lemma 4 that 2 a an d 21 are rational fractions of A, therefore they belong 
to Q (A) and then Q (p) = Q (A). □ 


3. From the Geometrical Property to the Dirichlet Properties 

In this section, the Geometrical Theorem is admitted, and we prove the 
Theorem 2 on Dirichlet Properties. 

3.1. Proof of Part a) of Theorem 2. 


Lemma 6 (Prism Lemma). Let g ^ be the column vectors of the 

matrix G( s ) generated at the s-th stage by the Smallest Vector Algorithm. 
Let the sets be the convex hulls: 


H'W = CONV ( gg (,s) , g } s) , g2 (,s) , 


We shall omit the indices Let H be the prism: H := H' + B. Then, 
with the usual notation h" for the orthogonal projection of h on B: 

-For each non zero integer point h in H , ||h"|| > ||gg|| holds. 
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-For each integer point h in H which is not of the form no go, with no 
integer, ||h"|| > ||g"|| holds. 

-For each integer point h in H which is not of the form nogo + nigi, 
with no and n\ integers, ||h"|| > 11 go 11 holds. 

- This implies that, i/(ho, hi, I 12 ) is a free triplet of integer points in H, 


then max I 

i=0,1,2 


|h"l 


> max I 

»= 0 , 1,2 


I Si 


= ||g£| 


Proof. Let h be any non-null integer vector in H. Since det (G^) = ±1, 
there exist three relative integers no, ni, n 2 such that h =nogo+^igi+« 2 g 2 - 
Let h' be the orthogonal projection of h on P. We have h'=nogo + n igi + 
n 2 g 2 and also h' G H 7 , which is the convex hull of six points. But an easy 
geometrical study shows that, in the plane, if p is in the convex hulls of six 
points, it has to be in the convex hull of three among these six points. In 
this proof, ’’positive” will mean: ”> 0”, and not ”> 0”. 

Then h' is the positive barycenter of three points among g' 0 , gj, g(,, — g' 0 , 
—gj, — g' 2 . This means that h' is the positive barycenter of three points of 
the shape £ogo, Mg}, £ 2 g 2 , with e, G {—1,1}. Then there exist positive 
real numbers t/o, y\ , t / 2 such that yo + Hi + V 2 = 1 and such that 

IT = yoeog'o + yi^ig'i + y 2 £ 2 g' 2 = nog'o + ngj + n 2 g' 2 - 

Then (n 0 - y 0 £ o) go + (n i - y\£\) g, + (n 2 - 2 / 2 e 2 ) g 2 = 0. But we have 

also 11X|| ||bg( s )|| ggT + ||X|| ||b"( s )|| g^ + ||X|| 11b^| g^ = X, like in the 
first Lemma of the Section 2, and then, with Zi = ||X|| ||bg( s )||, we have: 
z 0 gf Q + zigf + a: 2 g 2 = 0, with Zi > 0, by orthogonal projection on P. 

By the uniqueness of the barycentrical coordinates, up to a multiplicative 
coefficient, we get: for some real number A, 

no ~ yo£o = A20; ni - y X £i = Xz x ; n 2 - y 2 £2 = A z 2 . 

Then no — yo£o, n x — y\£ Xi n 2 — 7 / 2^2 have the same sign (0 has both signs). 
Without loss of generality, this sign may be supposed to be positive. Then 
we have n 0 > yo^o > — 1 , m > y\£ X > — 1 , n 2 > 7 / 2^2 > — 1 . 

- If for every i = 1,2,3, we have n* > —1 h, then no, n x , n 2 have the 
same sign. 

- If now, for some i, we have n x = !){£{ = —1, then ip = 1, and y.j = 7/^ = 0, 
with {i,j, k} = {0,1, 2}. Moreover, A = 0, and then n t = —1, rij = yj£j = 
0, np = y k £k = 0. Then, in every case, no, n x , n 2 have the same sign. Then 
||h"j| = ||n 0 go + nig" + n 2 gg || = |n 0 | ||goll + |ni| ||g"|| + |n 2 | ||g 2 ||, 

and the conclusions of the theorem become obvious. □ 


Let s notice that the Prism Lemma shows that in some way, our algo- 

(s) 

rithm gives best integer approximations of the plane on P. In fact, each g 0 
is a best approximation, and so are, after the nogg S \ the vector g^ and, 
after the vectors nogo’'* + n ig]_ S \ the vector gj, s \ The vectors g[, s \ g[*\ go^ 
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are not necessarily successive best approximations with disks of the eu¬ 
clidean norm in P, but so are they for the ’’hexagon” (which can be a 

parallelogram) H 7 ^, which depends on (gg S \ ) themselves. 


Proof. Let’s now prove the assertion a) of Theorem 2. 

We recall that, by convention, g"[ ) s+1 ^ 

we denote by gn , gm) the permutation of (gjf\ g^, g^ 


< 


„/d s +i) 

O 1 


< 


6 2 


and 
such that 


f(s ) 

gi 


< 


/(s) 

gn 


< 


/(s) 

p* v 7 

&III 


We admit the Geometrical Theorem, namely that there exists an infinite 

= L < Too. This will be 


set S of natural integers such that: sup 

seS 


'Ml 
| Siii | 

TTT 


proved in the next section. 

For any s € S', we consider the cylinder with center at 0, with basis in P, 
radius p ' 1 , and height 8 - 2 , namely: 


rM = Disk ' (T) + Disk "( 4 ') 


(the second disk being in B, and being a segment). The volume of T^ s ) 
is 8. Then, by Minkowski’s first Theorem, there is an integer point h / 0 
in For this theorem, see for instance, [6](Cassels), Theorem IV, page 
154. Then we have, by the Prism Lemma for the second inequality, 



or 




< M a constant, and the main statement of a) is 


proved. 

namely 


This last resu 


lim 

s —sGS* 


/(s) 

pp v ' 

&III 


t, 


with the help of the Lemma 2 of Section 2, 
= Too, implies the last statement of a), id est 


lim 

S —k + OO, S£S 



= 0. 


□ 


3.2. Demonstration of assertions b) and c) in Theorem 2. We shall 
need some well known results in Diophantine Approximation or Geometry 
of Numbers. 


Lemma 7 (Transference Theorem in dim 2). Let 1, cc, /3 be three rationally 
independent real numbers, let X be: X = 1 (1 ,a,/3) and let B be B = 
MX. Let h : = ( m,n,p ) be an integral point in Z 3 \ {0}. The following six 
assertions are equivalent: 

(a) inf IX • hi x [max (|?r| , |»|)1 2 > 0 

h =(m,n,p) with (n,p)^( 0,0) 
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(b) inf I ml x [max (lam — n\ , \Bm — pi)] 2 > 0 

h =(m,n,p) with m^O 

( a') inf |X • h| x [max (| ??7.| , |n| , \p\)} 2 > 0 

h=(m,n,p)£Z 3 \{0} 

( b') inf max (|m| , Ini , |d)x [max (\am — n\ , \/3m — p\)] 2 > 0 

h =(m,n, p)gz 3 \{0} 

(A) inf IX • hi x llhll 2 > 0 

hGZ 3 \{0} 

(B) inf ||h|| x ((am — n) 2 + (Bm — p) 2 '] > 0 

h=(m,n,p)GZ 3 \{0} V / 

Proof. The equivalence of (a) with (b) is a classical result. See for instance: 
[6] (Cassels), Theorem II and Corollary. The other equivalences are also 
classical and easy. □ 

Lemma 8 (Transference Theorem, geometrical point of view). Let 1, a, /3 
be three rationally independent real numbers, let X be: X = r (l,a,/3) and 
let O be ID = MX. As usual, let h' and h" be the orthogonal projections 
of h on D and P = D -1- . The following three assertions are equivalent, and 
also are equivalent to assertions (A) and (B) of the previous Lemma. 

(C) inf ||h"|| x llhll 2 > 0 
hez 3 \{o} 

(D) inf llhll x ||h , || 2 > 0 
hez 3 \{o} 

(E) inf ||h"|| x IIITII 2 > 0 

hGZ 3 \{0} 

Again, the proofs are easy. 

Definition. When one of the assertions (a), (b), (a’), (b’), (A), (B), (C), 
(D), (E) of Lemmas 7 and 8 is true (id est, all of them), it will be said that 
the couple (P, ID) is badly approximable. 

We shall also need famous Minkowski’s theorem on successive minima. 

Lemma 9 ( Minkowski’s Successive Minima Theorem). For any convex 
set E in M 3 which is symmetric about 0, let Ai(E), for i = 0,1 or 2, 
be the lower bound of the numbers A such that XE contains {i + 1) lin¬ 
early independent integer vectors. Then, if the volume of E is 8, | < 
Ao (E) Ai (E) A 2 (E) < 1 holds. See Cassels [5] (Cassels) Ch. VIII, page 
201 and following, especially assertions {12} and {13} page 203, or [6j 
(Cassels), Theorem V page 156. 

Proof of the b ) and c) of Th. 2 (Dirichlet Properties). We suppose 
that the Hypothesis of the b) Property of Theorem 2 holds: there exists 
c > 0 such that for any non-null integer point h, ||h|| ||h ,, || > c. Then, by 
our Geometrical Transference Lemma, there exists d > 0 such that for any 
non-null integer point h, 11 h'11 2 ||h // || > d. 
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Like above, let Tr be the cylinder: Tr = Disk 7 (R) + Disk 77 . 

( 7vd \ 3 2 

— J . The inequality ||h 7 || ||h 77 || > d implies that 

there’s no non-null integer point in AgT r. Then, for any R, Ao (Tr) > Ko¬ 
hl addition, we have Ao (Tr) Ai (Tr) A 2 (Tr) < 1 and also Ai (Tr) > 
Ao (Tr) > A'o , so that we get KqA -2 (Tr) < 1 and then A 2 (Tr) < K- 2 , 
2 

with K 2 = — 77). Then, for each R, the cylinder I\- 2 T r contains a free triplet 
Kq 

of integer vectors (ho, hi, h 2 ). 

IL'MII p(s) 

Let now be s in S, verifying ^ ^ ^ < L. We may suppose R = so that 

( 4A" 3 \ 

K 2 Tr = Disk 7 ( p + Disk 77 I -], which contains a free triplet of 

V 71 " ( p (s) ) J 

integer vectors (ho, hi, I 12 ), but the basis of which, Disk 7 (p^) is contained 
in the ’’hexagon” H 7 ( s ) generated by our algorithm. Then, by the Prism 
Lemma, 

Siii S 2 (S) <L 2 (p {s) ) 2 g 2 (s) < L 2 (p(*)) 2 f max h" (s) ) < 

yz=U,l,z J 

\K^1? 

--— and the main statement of b) is proved. The second statement fol- 

7r 

f(s) 

lows from this very result and from lim g Il: j- = +00 from Lemma 

s —^-)-oo, s£S 

2 at the beginning of Section 2. 

Let’s now verify the c) Property. We’ve just proved that under the 
Geometrical Theorem, for some M, sup ( gjjj g " 2 d J A M holds. 

seS \ J 

Now, if e s = det (B^) = det (G^) = ±1, and if (■ i,j, k) is a direct circular 
permutation of (0,1,2), forgetting the indices, we have: 

hi =e (g j A g fc ); b 7 = e (g 77 A g' k + g 7 A g"); b 77 = e (g 7 A g 7 fc ). 

Then for each i, ||b 7 || < 2||gj II || ||g 77 2 1|, and ||b 77 || < ||gj n || 2 . Then 

C=o a i X 2 b * f ' S) ) C=o aX 2 b,,?W ) - 4 ( g iii Ils2 (a) ||) < 4M 2 , and 


the main conclusion of the Part c) is proved. Furthermore, we have seen 
that for each s and each i, b 4 ^ < 2 g^ g"^ -1 ; in addition: 

gjjj S -2 ^ — M holds. Then, for each seS, hf s ' > < 

s') 

But, by the Lemma 2 of §2.1, lim g Ii:[ ; = + 00 . 

s —>-+oo 


2 M 
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Then lim 

s—¥+oo,sGS 


( max 

b' (s) 

Vi=0,l,2 

i 


0 and the second part of c) is proved. 

□ 


4. Demonstration of the Geometrical Theorem of §1.3. 


The demonstration of the Geometrical Theorem, Theorem 3. in Subsec¬ 
tion 1.3., involves only very elementary geometry, but is a little long. In 
order to prove it, we first need some auxiliary sets and definitions. 


4.1. The area and the set T (advances of 


>(s) 

pr v 7 

&III 


)• 


Notation 2. * We'll denote by twice the area of the triangle 


r(s) l(s) l(s ) 

g 0 Si g2 


* We’ll denote by (si S \sii\siii) the permutation of (g^, g} s \ 


such that 


r(s) 

Si 


< 


l(s) 

Sii 


s such that 


/(s) 

pr v 7 

&III 


< 


f(s) 

p. v y 

&III 


*In the Smallest Vector Algorithm, let’s denote by T the set of all integers 


< 


'(*+!) 

gin 


* The set T is infinite, because by Lemma 2, we have: lim 

s—>-+oo 

f(s) A f(s) , r(s) A f(s) , f(s) A f(s) 

gi A gn + gn A gin + gin A gi 


f( s ) 

p. v y 

&III 


= +OC. 


Remark, a) ^4^ = 

b) It’s easy to establish that for some i and j among {0,1, 2}, we have 

c) Then the sequence (yl^' ) ) A . gN is increasing. 


A( s+1 ) = y4( s ) + 


f(s) . /(s) 

g; A g, 


Notation 3. In this paper, for two non null vectors a and b in M 3 , we 
shall consider the angle of these two vectors corresponding to the canonical 
euclidean norm, and the measure of this angle which belongs to ] — 7 r; 7 r]. 
This measure will be denoted either by ”A (a, b) ”, or, more simply, when 
no ambiguity can occur, by ”(a, b) ”. 

Lemma 10 (Geometry on T). In the Smallest Vector Algorithm, for any 
s € T, ||gj^|| > - gjjj. . Moreover, for i, j € {0,1,2} we have: — < 


sif 

y(g'Kg?' ) 


Proof. Let i.j G {0,1,2}, i / j. We have 


/(s+l) 

gin 


< 


r(s) r(s) 

g \ y — o’ v 7 

l O'] 


then 


/(s+l) 

nr 7 

&III 


< 


'(*) 


+ 


T d s ) 


1 — 2 cos 


e' (s) e' (s) 

O % 1 Gj 


/(») 


'(*) 


and 


Without loss of generality, we may suppose 


g 


/(*) 


> 


g 


>V) 


. Then, by 


the cos formula above, if 


^(gf’.gf 


< ^ would hold, we would have 
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/(s+l) 

2 

/(a) 

2 

/(a) 

cr v ' 

&III 

< 

s i 

< 

nr v ' 

&III 


Then 


Moreover 


4 g ' w , g ?> 


/(a) 

nr v ' 

&III 


Hence 


/(s) 

gll 


< 


> u 
/(s+1) 

nr ' 

&III 


and s wouldn’t be in T. 
holds. 


< 


/(a) 

Si 


/(a) 

Sii 


< 2 


/(a) 

gii 


> 


/( s) 

pr v ' 

&III 


□ 


4.2. The notion of ”Needling”. As A.J. Brentjes has already noted 
in [3], if we want our vectors to have good approximation qualities, their 
projections g' 0 , gj, g(> on P must avoid the needling , i.e. flattening phe¬ 
nomenon. 

We’re going to define and study this phenomenon, but first we need the 
following Lemma in elementary geometry. It is very easy and its proof will 
be omitted here. 


Lemma 11. Let a and b be two vectors of the plane, with 0 < ||b|| < ||a||. 
Let's suppose that there exist real numbers s > 0 and M > 0 such that: 
e < \Z (a, b)| < tt — e and > M. Let p be the radius of the greatest disk 

centered at 0 and included in the parallelogram (a, b, (— a) , (— b)). Then 
there exists a real number M 1 > 0, depending only on e and M, such that 
> ML 

11 a 11 

As an immediate corollary, we have: 


Lemma 12. (Needling Sequence of Parallelograms) 

Let (a^ s \ b( s \ (—a^) , (—b( s ))) be a sequence of parallelograms, with: 
0 < ||bM|| < ||a( s )||. If these parallelograms are needling, id est if 


lim 

s —>-+oo 



0, where p is defined like in the preceding Lemma, then 


l i m sin(a^),b^M x 
s— H-oo V / 


bM 


= 0 . 


Proof. If the conclusion were false, then, for some p > 0 and for any s in 
some infinite set U, we should have: sin ^a( s ),b( s )^ > p and ---^y- > p, 


( s ) 

and then, by the preceding Lemma, |^ (s) y > M’ for some M’ and for 

s € U. But this negates the hypothesis of our Lemma. Then the conclusion 
is true. □ 


This last Lemma leads to the more important Lemma, which describes 
the needling phenomenon on the set T. 
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Lemma 13. (Needling Triangles) The three following assumptions are log- 

a( s ) 

ically equivalent: a ) lim -— — 2 = 0; 


s—>-+oo,sET 


■,(«) 


/(s) 

gill 


b ) lim -ji——77- = 0, where p (s ’ is the radius of the greatest disk cen- 


s—>-+ oo,sET 


f{s) 

p. v / 

&III 


tered at 0 and included in the convex hull: 
H /(s) = CONV (g^, gi, g' 2 , ~g'o, ~g'i, -g' 2 ) in P; 


c ) lim 

s—>+ oo,s£T 


UP 

Si 


/(s) 

Sn 


+ 7r — 


/(g£>,g S' 


= 0. 


Proof. First, let’s prove a)=^b). Omitting the indices I s !, we have: 

H' = cONv(gf ) ,g' 1 ,g£,-g' 0 ,-g' 1 - g^) inP. 

Let now K 7 be the convex hull of all the points 2g( — g 7 , with i ^ j and 
{i, j} C {0,1, 2} . Each g' belongs to K 7 , because 

s' = | (2g- g') + 5 (2g' - g') • 

In addition, the projection of the cofactors relation on P leads to: 

IIXII ||b"|| • g' + 11X|| ||b"|| • g' + 11X11 Ill'll • g' 2 = 0. 

Then, 0 is in the triangle ggg^g^- Let Fq be the homothety with center 
gQ and with scaling 2. 

Then Fq (0) = —gg is inside F (gggjg^), which is the triangle with sum¬ 
mits g 2 ; 2g , 1 — g^; 2g 2 — g' 0 . Then (—gg), and, in the same way, (—gj) and 
(—g 2 ), belong to K 7 . Then H 7 C K 7 ; then 7r (p^)~ < area(K'I s I). But 
K'(*> is formed with 13 triangles, each of them isometric to the triangle 
gggjgf>. Then: 7r (pD) 2 < area (K'^l) < 13ylI s I. 

Then, using a), we have 7r (p^) 2 < 

0. This implies b). Second, let’s prove b)=>c 


/(s) 

or v 7 
&III 


e(s), with lim e (s) = 

ir—>-+oo ,x£T 


If b) is true, we also have lim n p / ( (s > ,|i = 0, if p^ is the radius of 


s—>-+ oo,sET 


the greatest disk centered at 0 and included in the parallelogram with 
summits gjj, gj n , —g(j, — gj n . Then, by the last Lemma: 


lim 


f(s) f(s) . 

sm I gjj \ Siii ) x 


/(s) 

Sll 


/(s) 

Sill 


s— >-+oo,sGT 

T” of the last subsection above, 


= 0. But, by the Lemma ’’Geometry on 


g^ 


/(«) 

gin 


> i holds. 


— 2 


Then lim sin ( gn^,giii ) = 0. By the same Lemma: 

s—H-oo,s€T ' ' 


^(gp.gS 


> f. Then: lim 17r — 

,y s— >-+oo,s€T 


/(gp.gS 


= 0. 
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This last result, with the help of 


Lemma, leads to: limsup 

s— >-+oo,sGT 


4 (gf’.g?;’) 


4(g; (,) ,gS) 


> 4, by the same 


< We also have by the 


same Lemma, 


4 (gfU? 


> 5. Then: liminf 

' > s—>-+ oo,sGT 

n( s ) 


/(s) /(s) 

sm I g T Lgii 


> 0. 


But, if b) is true, lim n p /(s) ii = 0 holds, if p^ is the radius of 


s—>+oo,s£T llgjj^ 


the greatest disk centered at 0 and included in the other parallelogram 
with summits g{, g( n , -gj, -g^, and then, by 


lim 


pis) 




/(s) 

Siii 


> 4, we also have: 


„ 'I , = 0. Then, by the last Lemma on the needling parallelo- 

s—H-oo,s6T llgj/ ’ 


i- • / /(s) r(s) 

grams: hm sm g T , g TT 

s— »+oo,sGT ' 


we obtain: lim 

s— >-+oo,sGT 


gi (s) 


Si 


/(s) 


S '(s) 

&II 


= 0, and using the liminf above, 


& 


= 0. The proof of the part b)=>c) is done. 
Finally, the implication c)=>a) is obvious. □ 


In order to refute the needling phenomenon, we’re going to study what 
happens when the projections on P are ’’almost flat”. This will allow us to 
prove that the needling CANNOT happen. 


4.3. Almost Flat Triangles. Set T* of indices. 

Notation 4. T* will denote the set of all integers s € T such that the 


triangle gn^gi^g^* 5 ’’almost flat”, namely such that: 


4 (gS.gff) 


> 


30vr 

Hi' 


gi (s) 


1(b) 

Sii 


<0.1 and 


Lemma 14 (Flat Triangle Lemma). If s € T* (i.e. if the triangle 
go (s) g'i (s) g 2 (s) ’’almost flat”), then we also have the following relations: 


1571 

~3r 


< 


4igS.g ?*> 


177r 157T 

<- and - < 

31 31 


/(gSfg'P 


< 


and also: ■ 

f(s) 

Sii 

f(s) 

■ > 0.979, 

'(s) , '(«) 

&II ^ &III 

< 0.23 

/(s) 

pr v ; 

&III 


Shi 






177r 

Hi’ 


and: 


i(s) i(s) r(s) 

g \ / — nr v ' — nr v ; 

I &II &III 


< 0.33 


l(s) 
nr v ' 

&III 


Proof. Because in T, 
lowing inequalities also hold: 


'(s) _ Us) 

Oi o? 


> 


/(s+1) 

pr v ' 

&III 


> 


r(s) 
nr v ' 

&III 


/(s) /(s) 

pr v / — pj. v j 

&III &I 


> 


/(s) 

p. v / 

&III 


and 


/(s) l(s) 

Six - Si 


> 


r(s) 

pr v / 

&III 


> 


holds, both fol- 


r(s) 

Sii 
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The first inequality leads to: 


r(s) 

gin 


< 


l(s) 

gin 


+ 


l(s) 

Si 


/(s) 

&III 


g?*> 


/(s) 

gin 


- 2 cos (g'n^gi^ 


Then: 0.1 — 2cos ( gj-jj, ) > 0; then: cos ( , gj iJ ) < 0.05. 


/(«) _'(s) 


^(gS.gf* 1 


Hence 

n the same way, we obtain: 

and 


> -yph Hence 




< 


1771 


— 31 


^ (gS'.gf 1 ) 


> 


157r 


— 31 


4gS.gf*> 


Again, let’s use the inequality: 


/(«) 

o' v ' 

gin 


< 


l(s) 

gn 


+ 


g'l w 


/(s) 

or v ' 

&III 

2 

< 

2 cos | 

/(s) 

p. v / 

Oil ’ 


< 


17tt 


— 31 • 

/( s) r(s) 

gn “gi 


/(s) 

gn 


. Then 


r(s) 

Si 


But -2cos ((gi ( i ) ,gi (s) ) ) < -2cos (4f) < 0.31. 


Then, dividing by 


/(«) 


gn 


/(*)! 


1 < 


Sii 


e 'M 

sill 


4 + 0.01 + 0.031 


, we obtain 

, i.e.: x 2 + 0.031x — 0.99 > 0, 


+) 

Sn 


/(«) 

Sin 


with x = 


'M 

Sn 


-'<«) 

&III 


. This implies x = 


/(s) /(s) 

pr ' -L p 1 ' ' 

&II ' &III 


g'n } + gin 


< 


/(s) 

2 

/(s) 

gn 

+ 

nr v ' 

&III 


& 


'M 

&III 


> 0.979. We have 


+ 2 cos (j§r) 


'(+ '(«) 
gn ' gin 


< 2 

< V0.0521 


/(s) 

gin 


- 1.9897 x 0.979 


/(s) 

gin 


/(+ 

or v ' 

&III 


/(s) 

gn 


then 


, then 


< 0.23 

Hence the last conclusion is obtained. 


/( s ) 

pp v ’ 

gin 


/( s ) 

pp v y 

gin 


□ 


4.4. From a Lemma to the Geometrical Theorem (Theorem 3). 
Notation 5. The sequence (<V S )) +• defined by: <V S ) = T 1 . 

| g m | 

//limsup <V S ) > o, then the algorithm is balanced , i.e. the triangles on P 

s—>+oo 

do not needle. 

Lemma 15. (Monotonic Subsequence Lemma) Let [m;+oo[ be an 
interval of N such that [m; +oo[nT C T*, which means that every advancing 
triangle with its range in [m; +oo[ is almost flat. Then the sequence (aV^) 
is increasing on [m;+oo[nT, which means that for any s,t € [m;+oo[nT 
with s < t, j holds. 

t [S J M t) 

Siii Sin 



26 


Christian Drouin 


Let’s admit this Lemma, which is proved in the next Subsection. Then 
we can demonstrate the Geometrical Theorem, Theorem 3. 


Proof of Theorem 3. Let’s suppose that the conclusion of the geometrical 


( s ) 

Lemma is FALSE, namely that lim n p . .m = 0. Then, by the Lemma 

s — }-\-oo,sGT gm 


on needling triangles of the second subsection above, lim 




s—H-oo,seT|| g 'W 


= 0. 


By the same Lemma, we have 
+ (n - 


lim 

s— >-+oo,sGT 


/(s) 

Si 


/(s) 

Sii 


/(gS'.gs 


= 0. Then, there exists an in¬ 


teger m such that [m; +oo[nT C T*, which means that with a range great 
enough, any advancing triangle is almost flat. Then, by the Monotonic Sub¬ 
sequence Lemma, the sequence of the a^ ^ with s > m is increas- 


I[£ ,(S) I 

Bill 


A<» 


inq on T, which is infinite. This is contradictory with lim , = 0. 

s_ H- 00 > seT ||gin || 

> 0 and the Geometrical Theorem is 

□ 


Then, we have lim sup n p !,!u 

s— H-oo,seT Sm 


proved. 


This Monotonic Subsequence Lemma has now to be proved. 


4.5. Proof of the Monotonic Subsequence Lemma. Let s be an in¬ 
teger in the interval [m;+oo[n T C T* as in the Hypothesis. Let’s denote 
s' the successor of s in T, i.e. the smallest integer t in T such that s < t 
In order to establish our Lemma, it suffices to show that 


(4.1) 


A (*) 


At*') 


l(s) 

2 — 

/(s') 

2 

Bill 


&III 



We have four cases: 

Either g'^ 1 ’ = gi (s) - gii^ (Case (I 


/(s+l) /(s) 

or gill = gii 


l(s) 

Si 


\ 

(II x 


TT n /(s+1) f(s) Its) 

II)), or g^ = gi - gni 


/(s+l) /(s) 

or g m = g m 


(I N III) 

(III X I). 

We’re going to prove the inequality Hnea 14.ll) only in the case (II 
The demonstration is similar, or easier, in the three other cases. 
Case (II \ I): gm +1) = gf? - gj a) . 

We have to obtain first 


/(s) 

gi 


X I) 




Af( s+1 ) 




/(s) 

Sii 


/(s) 

gi 


/(«) 

2 — 

'(s+l) 

2- — A ( S ) 

'(s) 

2 

Bill 


&III 


&III 
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It’s sufficient to show: > | 1 + 


But 


_4(s+l) _ - 4(s) +||gm ) Agi (s) || 


g( W 


&III 


AM 

sufficient to show: 


AM 


= 1 + 


Ilg^Ag^l 


g ; (s) A g ; ( : s) + g g s) Agjw+ g &> Agi w 


It’s 


1 + 


||g^Agi (S) | 


|gi (s) A g n S) +gn S ^ Agj^+gjjj 5 A g j (s) 


fgglV'l 


11 gi (s) A g ii S) +gii S) Ag^ + g m Ag' (3) 11 - 


gi (s) Agj( s) 

• -1- ■ 

g^Agg? 

g^Ag( W 


g^Ag( (s) 


1 + 


enough to show: 1 + 



i.e.: 


+ 2 |; i.e.: 

. To obtain that, it’s 


gj (s) Agjj S l 

■ -U . 

gg'Wg? 


Sm || 

gi^Ag?* 5 

1 

gi^Ag?^ 

“ 2.1 

Nil 


i.e.: 


(4.2) 


2.1 

/(.s) 

Si 

2.1 

. L 

/( s ) 

gl 


f(s) A /(s) 

gi Agn 

2.1 

. L 

/(s) 

gl 


/( s) . f(s) 

Si! A gin 


/(s) 

pr v ' 

&III 

I 

/(s) 

nr v ' 

&III 


/( s ) . r(s) 

Siii A gj 


/(s) 

p. V 7 

&III 


f(s) A /(s) 

gin A gi 



< 1 


We 


nave: 


2.1 

l(s) 

gl 


/( S ) A /(s) 
gi AgjV 

2.1 

/(a) 

gl 

sin ^ 

(«?'>.«S>) 

) 

t(s) 

Sii 


/(a) 

nr v ' 

&III 


/( s ) A /(s) 
gin A gj> 


/(«) 

or 

&III 

sin ^ 

( g ; w ,gS) 

) 

e /(s) 

gin 


2.1 x 0.1 

< - TvFTL < °- 22 

sin (Vf) 

In addition: 


2.1 

/(a) 

gl 


/(s) 

or v ' 

&III 



2.1 sin (^zi) 

<- /1 V 7 3 . 1 ’ < 0.22. 

- sin(fi) - 


f(s) A f(s) 

gii Agm 


f(s) A f(s) 

gin A Si 


2.1 sin 


/(a) /(a) 

gm.gii 


sm 


/( s) l(s) 

Siii, Si 


Finally: 


2.1 X 

/(s) 

gl 


/(a) 

or v ' 

&III 



< 2.1 x 0.1 < 0.21. 


/( s ) 

gn 


/(s) 

nr v ' 

&III 


The three last inequalities lead to the sufficient condition: ilnea l4~2l) . 

4( s +!) 

'(s + l) II 2 ■ 


Then we have proved: M 3 ' s) j < j4(s+1) 
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If (s + 1) € T, the proof of 

< 


then we have both 

^ 4 ( 3 + 2 ) 

0>+57F 

&III 


/(s+1) 
qc v ' 

&III 


Inea 14.ip is finished. If now (s + 1) ^ T, 

l(s+2) 

cr v ' 

&III 


, and Al^O < then ^ ts) 

’ — '(a) 


Si 11 


72 < 


and again the proof of (Ineq ITTIl is finished. 


In the same way, as long as (s + i) (f T , for i = 1,2,..., we have: 


AW 


I '(*) 

Sin 


~2 < 


A( s + 2 ) 

I '(*+2)1 

&III 


7 < 


4(*+3) 


I '(*+3) 

SHI 


77 < ••• < 


^( s +*) 


II '(*+i) 
Bill 


77 < < 


A( s ') / u ■ 

TiTTF’ s bemg 


IIE ,(S,) 

&III 


the successor of s in T. 

AW 

Then 


l(s) 

2 — 

/(s') 

2 

&III 


Bill 



and the conclusion (Ineq lTTT) is reached in the 


case (II \ I). 

The reasoning is similar in all the four cases. 

Then the conclusion of the Monotonic Sequence Lemma is es¬ 
tablished. So is the Geometrical Theorem, and also the Dirichlet 
Theorem and the Lagrange Theorem, but the last one only in a special 
case. We have to prove it generally. 


5. Lagrange Theorem from Dirichlet properties: complete 

demonstration 

Now, using the Theorem on Dirichlet Properties, we prove the Lagrange 
Theorem with the help of some Definitions, Lemma, Propositions. First we 
give the statements, then the proofs. 


5.1. Definition and Statements of §5. 


Definition (Max-Dirichlet Property). It will be said that a sequence 

( pW ) = (pS” , p ( j' , ' ) , P 2 ^ of triplets of integer vectors has the max-Dirichlet 

Property concerning D = MX (resp: P = X 2 -) if there exists an infinite 

2 


subset 5 of N such that: sup 

s£S 

with 


max 

t=0,l,2 


As) 


max 

*= 0 , 1,2 


P 


'00 


lim 

max 

/(s) 

Pj 

s — >’H-oo,sE.S' 

U=0,l,2 



= 0 (resp: lim 


[ max 

p't> 

V 2=0,1 ,2 



< +00 , 

= 0 ). 


Lemma 16 (Polarity and Dirichlet Property). Let (P^) be a sequence of 
integer matrices, all with the same determinant D > 0, up to the sign, id 
est, for each s € N, det (p( s )) = e^D, with € {—1; 1}. Let (p( s ))* be 
the polar matrix o/P( s ). Let X be a triplet of rationally independent real 
numbers. 
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If the sequence has the max-Dirichlet Property for the plane P = 


X x then the sequence yD. (p( s )) 
line D = MX. 


has the max-Dirichlet Property for the 


Proposition 1. Let 0 be a real root of a third, degree irreducible polyno¬ 
mial P (t) = t 3 — mt — n, with m and n rationals. Let the vector © be 
© = T (1 , 9, 6 2 ). Let R be any rational matrix, with det (R) f= 0, such that 
CR has integral coefficients, with C an integer. Let X be: X = R© = 
R T ( 1,9,0 2 ). Let’s suppose that: X = 1 (xq,xi,X 2 ), with 0 < xq < x\ < 
X 2 ■ Let be the sequence of integer matrices generated by the Smallest 

Vector Algorithm with initial value X = R0. 

Then the sequence (G^)) has the max-Dirichlet Property for the approx¬ 
imation of P =X J - and ^C. T RG^^ has the max-Dirichlet Property for 

the approximation of II = ©^. Moreover, there exists an integer A such 
that the matrices A (s ) = AR 1 (G( s ))* are integer, and such that the se¬ 
quence (A (s )) has also the max-Dirichlet Property for the approximation of 
A = M0. 


Proposition 2. Let 9 and © be like in the previous Proposition. Let A be: 
A = M@. Let (A( s )) be any sequence of integer matrices having the max- 
Dirichlet Property for A, and all having the same determinant D > 0, up 
to the sign. Let (J^) be the polar matrices of the (A^). Then there exists 
a sequence of rational matrices and an integer Q, which depends 

only on m and n, such that: 

• For each s, © is an eigenvector for M( s ); 

• lim inf (||Z)Q 1 M (s )j( s )||) < +oo, the matrices (DQ 1 M^j( s )) hav¬ 
ing integral coefficients. 


Proposition 3. Let 9 be, like in the two previous propositions, a real root 
of a third degree irreducible polynomial P (t) = t 3 — mt — n, with m and 
n rationals. Let X = r ( xo,xi,X 2 ) be a free triplet of three positive real 
numbers from the ring Q [9\. Then the Smallest Vector Algorithm applied 
on X makes a loop: there exist integers s and t, s / t, and a real number 


A such that: 


( 

_//(«) 


„//(») 


g 0 

5 

g 1 



( 

a-" W 


„//(*) 

\ 

g 0 


g 1 



5.2. Demonstration of the Lemma. Let’s denote: 

(s) (s) (s) 

qo \qi ,02 


(p(*))* = q(») = 


; then, for any direct circular permu¬ 


tation ( i,j,k ) of (0,1,2), forgetting the indices W we have: 


q* = % (Pj A p fc ); q i = %( P j A p), + p' A p") ; q" = fJ p' A p[ 


Let’s denote: max 

*= 0 , 1,2 


P 


/(s) 


= max 

*= 0 , 1,2 


P 


i(s) 


= 11 "W. 
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By hypothesis, (//( s )) 2 < L holds for some L and for every s in the 


infinite set S. For each i, and any s G S, we have: 


and 


/(s) 


max 
*=o,i, 2 


<1 


'(«) 


< (^ /(s) ) 2 - Then 


l(s) 

q> 


d q : 


/(s) 


max 
*= 0 , 1,2 
2 


q'f’ 


< 


< 2. .As) Ms) 


max 

2 — 0 , 1,2 


Dq'T 


^((a/( s ))V ( s) ) < -pg-, and 

< 4L 2 . 


then: max 

\*=0,i,2 

We have to prove in addition that the limit of the first factor is null. 


lim 

max 

p"S s) 

s —H-oOjSES 1 

l 2=0,1,2 



lim 

max 

/(s) 

p, 

s — 

U=0,l,2 



= 0. This implies: 

= Too. Otherwise, the set of all the integer 


(s) • • 

vectors p. , with s in some infinite set T C S, would be bounded, and 


( max 


V 2=0, 1,2 



would have a non-null 


sGT 


minimum. Contradiction! 

Then we have: lim fi'^ = Too. But we also have, for every s € S, 


s — y~\~oO)S£S 


q, 


'(») 


< - 2 . V 

—h ~ d 


n{s) 


Then lim ( max 
s—>-+oo,sES' V 2=0,1,2 


D q' 


T 

/(s) 


/(«) 


< 


2 L 


1 


D ‘ 


= 0. 


We have established that the sequence yD. (p( s ))*j has the max Dirich- 
let Property for the line B = MX. 


5.3. Demonstration of Proposition 1. Let 0 be a real root of a third 
degree irreducible polynomial P (t) = t 3 — mt — n, where m and n are 
rationals. 

For the initial value X = R0 = R 1 (l, 8, # 2 ), let (B^) and (G( s )) be the 
sequences of integral matrices generated by the Smallest Vector Algorithm. 

First, we establish that the couple (P, B) = (X 1 ,KX) is badly approx¬ 
imate, in the sense of the Lemma 8 and the following Definition in Sub¬ 
section 3.2. 

By a classical theorem that we have already cited, (see [6] (Cassels), 
Theorem HI, page 79, statement (2)) the couple (@- L ,M0) is badly ap¬ 
proximate. Then: inf |k*0| 

k integer /0 L 

Let’s suppose that the couple (P, B' 
imable. Then we would have: inf 

h integer 7^0 


> 0 . 


= (X i ,KX) is NOT badly approx- 


|h • R0| . Ilhll 


= 0 ; 
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then inf I ( 1 Rh) • ©I . II 1 Rh||~ = 0. Let q be an integer num- 

h integer jtO L J 

ber such that qH has integral coefficient; then 

inf | (q 1 Rh) • ©| . ||« 1 Rh||“ = 0, with (q 1 Rh) non-null inte- 

h integer ^0 L 1 V ' 111 11 J V ' 

i 11 112 

ger. Then inf |k • ©| . ||k|| =0. Contradiction. Then the couple 

k integer ^0 L J 

(P,B) = (X^MX) is badly approximate. 

By the Dirichlet Properties Theorem of Subsection 1.3., part b), the 

sequence (Gl s )), generated from X by the Smallest Vector Algorithm, have 

the max-Dirichlet Property for the approximation of P =X~ L . There exists 

an infinite subset S of N such that 

r / . ,\/ ..,\ 2 i 


sup max g,; 
seS V=°V 2 


max g. 

*= 0 , 1,2 l 


< +oo, 


with lim I max g, 

s—>-+oo,sES' V i=0,l,2 


sup max ( C. 1 Rg - 

sSS y*=0,l,2 I V 


with lim f max [C. x Rg- 

s-H-oo,seS V *=0,1,2 V ' 


R0 =0. Then: 


max \\C. 1 Rg 
*=0,1,2 II 1 

)•&{) = 0 . 


< +oo, 


That means that the sequence yO. 1 RG^j has the max-Dirichlet Prop¬ 
erty for the approximation of II = 0^. 

But the (yO. r RG^^ have all the same determinant (G 3 det (R)) , up to 
the sign. Let A be A = C 3 det (R). Then, by the previous Lemma, the 
sequence ^A^RG^^ ^ = ^AR -1 (G^)*^ of integer matrices has the 
max-Dirichlet Property for the approximation of A = M0. 


5.4. Demonstration of Proposition 2. Let (A^) be a sequence of 
integer matrices having the max-Dirichlet Property for A = M© = 

M T (l, 9, 6 2 ), with # 3 = md + n. 

We suppose that the matrices A^ have all the same determinant D > 0, 
up to the sign, which means that for each s € N, det (A^) = s^D, with 
€ {— 1 ; 1 } . 

Let ^a^, a[ s \ aj, 5 ^ be the column vectors of A^ s \ We choose one of 

these three vectors, say ag A \ which will be more simply denoted: a^ := 

a ( j S ' 1 . Let’s define its coordinates by: a( s ) = 1 cty S \ . 

Let’s denote: ii' 1 " 3 ' 1 = max a^ A ^ and = max a"j A ^ . We sup- 

*= 0 , 1,2 1 *= 0 , 1,2 * 

pose that there exist an infinite set S of integers and a real number L such 
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that for each s € S, the inequality n"^ < L holds. Let s be any 

element of S. From now on, we may omit the indices ( s ). 

The notation M = M 1 " 5 ) will denote the rational matrix 
( Ttl.dx T &Z v dx \ 


= I n ■ a x a z a y | . it is a natural sucli tnat 

y n ■ a y n.a x + m.a y a z ) 

Qm and Qn are integers, then has integral coefficients. 

We have: M© = M ( 9 ] = (—m.a x + a z + a v 0 + a x 6 2 ) ( 9 | . 

U / " U 7 

Let’s denote by A the following element of Z [9]: 

A := (— m.a x + a z + a y 0 + a x 9 2 ). Then we have M© = A©; A is an 
eigenvalue of M with eigenvector © . 

Let the polar matrices of the (A^). We consider the sequence of 

the matrices 11^ = II = i MJ. 

Let (a##, a#, a) be the three column vectors of M. Then, with scalar 

/ a ## *j 0 a ## *ji a ## «j 2 \ 

products: 11= 1 MJ= a # «j 0 a#«ji a # *j 2 = say, 

\ a«jo a • ji a • ji / 

with i = 1,2, 3; j = 1, 2,3. 

We now have to find an upper bound for each of the 

l a r \ l 0 1 0 \ 


. If Q is a natural such that 


products: II = 1 MJ = 


= 0 ij), say, 


We have: a# = Q# 


,## _ 


= Q## 


= Q#a, with Q# = 


= Q ## a, with Q** = 


/ 0 1 0 
= 001 
y n m 0 
—m 0 1 \ 

n 0 0 J 

0 n 0 / 


We have: Q#0 = 00, and Q##0 = (-m + 6 2 ) ©. 
First let’s consider the la# • i,|. 


Let v be: v : = - 


Then also Q#V =9v. 


With always the same kind of notations, we have a# = a#' + a#", and: 
a#•jj = (a #/ + a#")«(j' + j") = a #, «j•+a#"«j' / = a #/ «j•+a#"«j", with 

j % = % ( a j A a fc ); j - = £ fa" A a' fc + a' A a") ; j " = % fa'- A a' ) . Then: 


(5.1) |a#.j, t |<l(||a#'|| n"n' + + ||a#'|| /iV / + ||a#"|| ( M ') 2 ) 

and we have also a# = Q#a = Q# (a^ + ||a^|| v) = Q#a(, + ||a^|| 9u. 
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This proves first that the distance ||a 7 ^ / || between a 7 ^ and D is less than 

|Q #a o|| : 


(5.2) 




< 


Q #a o 


< 


Q 


# 


X llaj < 


Q 


# 




(this using the norm of the matrix). 

Moreover, we have a 7 ^' + a#" = a 7 ^ = Q#a() + 11aQ 11 0v 
then: a#" = 11a^11 Qv + Q^ap— a#'. Then: 


(5.3) 




< K 


+ 2 


Q 


# 


X 11 a 0 11 < /r + 2 


Q 


# 




Putting 15.21 and 15.31 in 15.11 we obtain: 


a #( s ) • j) s) 


and then: 


< b iv ,{s) ) F" {s) (2 ||Q # || +0) + i ||Q # || {fi' is) y 




L , 

s id 2 


Q 


# 




Q 


# 


A* 




The limit of the last term is 0. 


Then the set of the 


»#(«) • 


, with s in 5, is bounded. A similar 


demonstration shows that the 


(s) .(s) 

*o «Ji 


,##(«) . i( s ) 


are also bounded and so are 


Then the set of the LE^ = 1 M( s )j( s ) i 


is 


in an obvious way the 

bounded. But = (A^)*, with det (A^) = ±D. Then the matrices 
have integral coefficients. We have seen that the matrices 
have also integral coefficients. In addition, the sequence (DQ 
is bounded, and the proof is done. 


5.5. Demonstration of Proposition 3. Let 6 be a real root of a third 
degree irreducible polynomial P (t) = t 3 — mt — n, with m and n rationals. 
Let X = T (xo, xi , X 2 ) be a free triple of three positive real numbers from 
the ring Q[0], Let 0 be © = 1 (1,0,0 2 ). Then there exists a rational 
matrix R , with det (R) 7 ^ 0, such that X = R0. 

Then, by Proposition 1, there exists an integer A such that the matrices 
A^ = AR _1 (G^)* are integer, and such that the sequence (A^) has the 
max-Dirichlet Property for the approximations of A = M0. All the matri¬ 
ces AW have the same determinant, say D > 0, up to the sign; then, by 
Proposition 2, there exists a sequence of integer matrices (M( s )) such that, 

for each s, © is an eigenvector for and lim inf (IlDQ^M^J^II) < 

s— H-oo 11 J 

+ 00 , the matrices (DQ. T Ml s l having integral coefficients, with 

j(A = (AW)* = (AR - 1 = A -1 t RG ( s ). 

There exists an infinite subset S of N, such that the set of all the integer 
matrices ( DQ with s in S is bounded; then it is finite. Then 
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there exist s and t, t / s, such that = r M^ T R G^K By 

transposition: (B^) 1 RM^ = (B^) 1 RM^. We apply that to the 
column vector ©: (B^) 1 RM' S @ = (B^) 1 RM^0; then, by ’’eigen¬ 
vector”, and because R© = X, we have A (s > (B^) 1 X = A^ (B^) 1 X, 

then (B^) _ 1 X=^(bW)" 1 X. 

This reads: X* s ) = AX^\ with A = and our Lagrange Theorem is 
proved if X = T ( xq,x\,X 2 ) is a free triplet of three positive real numbers 
from the ring Q [ 0 ], 0 being a real root of a third degree irreducible poly¬ 
nomial P (t) = f 3 — mt — n, with m and n rationals. Of course this case is 
general, as we’re going to verify it. 

5.6. Prom Proposition 3 to the Lagrange Theorem. This part is 
very quick. Let p be a real root of a third degree irreducible polynomial 
S ( t ) = f 3 — at 2 — bt — c, with with a, b, c rationals. Then 0 = p — | is 
a real root of a third degree irreducible polynomial P ( t ) = f 3 — mt — n , 
with m and n rationals. If xo, x±,X 2 are elements of Q [p], they also belong 
to Q [0] = Q [p]. Then Proposition 3 implies the conclusion of the main 
Lagrange Theorem (first part). The second part of the theorem has been 
established in Section 2. 

6 . Bibliography and Themes related to this Paper 

The work nearest to the present paper is the book by A.J. Brentjes [3j. A 
lot of themes are in common: the approach of the continued fractions with 
matrices and linear algebra, the fact that non vectorial algorithms are used, 
the study of angular properties and of the needling phenomenon...Brentjes’ 
book is mainly concerned with algebraic results, best approximation, and 
(strong) convergence, rather than with ’’Dirichlet” approximation, with the 
optimal exponent, or ’’Lagrange” results. However, it contains a Lagrange- 
type statement, in the Corollary, page 106, but with a lattice which is not 
Z 3 . 

Most multidimensional continued fractions algorithms, among those 
which are additive (or subtractive, or multiplicative), are of the vectorial 
type. This means that in such an algorithm, the vector Xl s+1 ) depends 
only on X( s ) = ||X|| ^gQ^, g^ s \ j i n a simple way, and not on 

^gQ S \ g^ b \g'j in the plane P. In this case, the algorithm defines clearly 

a discrete dynamical system, the orbits of which are the sequences (X( s )). 
There are a lot of interesting studies of these dynamical systems, by Fritz 
Schweiger, J.C. Lagarias and many others, but our algorithm, like Brentjes’ 
one, is non vectorial, and different techniques are used. With such non 
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vectorial algorithms, results everywhere may be obtained. With vectorial 
ones, most of the results are obtained almost everywhere. 

Apart from Brentjes’ book, there is another treatise by Fritz Schweiger on 
multidimensional continued fractions [22]. It is very complete and presents 
the general Brentjes’ algorithms, but deals mainly with vectorial algorithms 
and dynamical systems. 

The continued fractions are only a tool in the theory of Diophantine 
Approximation. Here we use two theorems by Minkowski in Geometry of 
Numbers. The references in these fields are for instance: 0,0, m , and 

m- 

In the area of algorithms which aim to best approximation, apart from the 
specific Brentjes’ algorithm, we may cite the Furtwangler’s algorithm [12] 
(an error was pointed out by K.M. Briggs, see his paper), which inspired 
Keith Briggs 0] and Vaughan Clarkson [91; see also the Ph. D. thesis of 
V. Clarkson: 0. 

There are some studies of the matrices of best approximations, which 
could be connected to our work: By J.C. Lagarias: [18] . and and a 
review by N.G. Moshchevitin: [20] . 

In the present paper, the result on best approximations is the Prism 
Lemma, at the beginning of Section 3. It is an easy result, but perhaps 
it clarifies the problem. It is more efficient if the hexagon it involves is 
balanced, and we have some results in this direction in this paper. 

J.C. Lagarias has also build in [19] a very interesting algorithm, which 
is additive but not positive, and which provides best approximations. See 
also the very complete paper by N. Chevallier: [7]. 

The LLL algorithm (named after A.K. Lenstra, H.W. Lenstra, L. Lovasz) 
is very efficient in Number Theory. It provides good approximations, and 
even Dirichlet approximations, with the optimal exponent: see [2], by W. 
Bosma and I. Smeets. But maybe it is not designed to possess approxima¬ 
tion properties with triplets of integer vectors, nor Lagrange properties, as 
the Smallest Vector Algorithm does. 

There is an another kind of Multidimensional Continued Fractions, very 
different from the additive (i.e. subtractive) ones we have considered un¬ 
til now. These other constructions use stars of sails, obtained from hy¬ 
perplanes and pyramids in M fc . The original idea is due to K. Klein, H. 
Minkowski, and G. F. Voronoi. V. I. Arnold renewed the interest toward 
this theory: [1], ’’Lagrange” results seem to have been obtained, by G. 
Lachaud, [16], E. Korkina m, or O.N. German and E.L Lakshtanov: m- 
But their statements don’t seem as simple as the Theorem 1 of the present 
work. In the cited paper, V.I. Arnold has written: 

’’The attempts to generalize to higher dimensions the algorithm (empha¬ 
sized by V.I. Arnold) of continued fractions lead to complicate and ugly 
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theories. For instance the sail corresponding to a cubical irrational number 
is a double-periodic surface. However the algorithms define instead of this 
surface a path on it. [...] the path is not periodic at all and looks like a 
rather chaotic object; it is unclear how to describe the cubic irrationals in 
terms of the combinatorics of this path”. 

We can just hope that Arnold was only partly right. It would be inter¬ 
esting to study the relation between the regularities we have pointed out 
in the ’’chaotic” paths generated by our algorithm for cubic numbers, and 
the symmetries of the corresponding sails. 
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